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■ The recently found shock wave solution in the scalar field model with 

\Q ' the field potential V(<p) — \<j>\ is generalized to the case V(<j>) = \4>\ — \\<t> 2 . 

f^*) , We find two kinds of the shock waves, which are analogous of compression 

CN ' and expansion waves. The dependence of the waves on the parameter A is 

investigated in detail. 

PACS numbers: 03.50 Kk, 05.45.-a, ll.10.Lm 
p 1 1. Introduction 

fH \ Interesting and still poorly studied group of field-theoretic models are 

these with V-shaped potentials. Physical systems with just few degrees of 
freedom and V-shaped potential are studied quite frequently. There are nu- 
merous results for systems such as e.g. bouncing oscillators. These systems 
are mainly studied in context of chaotical behaviour and grazing bifurcation 
[1-7]. The V-shaped potentials appear in research of plasma physics [8]. Fur- 
thermore, they can play an important role in pinning phenomena which can 
describe a process of vortices attaching to lines of impurities [9,10]. Apart 
from the applications, they are also very interesting on purely theoretical 
grounds because of scale invariance of vacuum sector, [11]. 

There are only few analytical results for the field-theoretic models with 
V-shaped potentials. An example of such model, which originally derives 
from a mechanical system, has been proposed in [12]. The model considered 
there is obtained as a continuum limit of the system of coupled pendulums 
that are allowed to take the angular positions belonging to the interval [— </>o, 
4>o]. This leads to appearance of a V-shaped potential, and consequently to 
nontrivial dynamics of the system. The angle </> = corresponds to the 
pendulum in upward vertical position. A simpler model has been proposed 
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in [11]. It derives from coupled system of bouncing balls. The potential in 
that model has the form V{<f>) = \<ft\. It can be regarded as a limit case of a 
large group of V-shaped symmetric potentials. 

In the present paper we consider the classical scalar field model with the 
potential V(4>) = \<ft\ — \\<ft 2 , where A is a real constant. Such potentials with 
A ^ naturally appear in the case of system of coupled pendulums [12], as 
well as in the system of bouncing coupled balls obtained from system studied 
in [11] by adding new couplings, see Fig. 2. Specifically, we investigate shock 
wave solutions^ In the particular case of A = they have been discussed in 
[13]. 

The field-theoretic model with the potential \(j>\ has special symmetry. 
If 4>(x,t) is solution of equation of motion, then u 2 c/)(^, -), where v is a 
positive constant, obeys this equation as well, see [13]. This symmetry is 
"on shell" type because the action functional is not invariant with respect to 
the scaling transformation. In general, in real physical system, apart from 
term |0|, the V-shaped potential has also another terms. A squared term is 
an example of the simplest perturbation that breaks the scaling symmetry. 
Investigation of effects of such perturbation is a very important issue. 

Our paper is organized as follows. In Sec. 2 we show connections between 
the model and physical systems. Next section contains analysis of solutions 
that have the properties of expansion shock waves. In Sec. 4 we investigate 
symmetric compression shock waves. Sec. 5 summarizes the paper. 

2. The modified model and related physical systems 

2.1. General definition 

Let us consider the model of a real scalar field (ft in one dimension. The 
Lagrangian has the form 

L=\{d t cftf ~\{d x cftf -V(cft), (1) 

where <p = (ft(x, t) depends on rescaled position x and time t which are 
dimensionless. The potential has the form 

v(<f>) = \<ft\ - \\4> 2 . (2) 

1 In this and previous papers we use terminology " shock wave " so as to name disconti- 
nuity that moves. There is a criterion in theory of hydrodynamics that distinguishes 
between shock waves and different kinds of discontinuities as e.g contact discontinu- 
ity, tangential discontinuity etc, see e.g. [15]. Our system is clearly different from 
e.g. gas medium and the criterion can not be applied directly. Every choice of name 
for discontinuity in model considered here can merely reflect some kind of analogy 
between the discontinuity we describe and hydrodynamical discontinuity that has 
this name. For this reason we shall remain at the name " shock wave " and discard 
looking for a better term. 
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The evolution equation which corresponds to Lagrangian (1) has the form 

(^ 2 -^)0 + sign(0)-A0 = O. (3) 

Because of sign function equation (3) is clearly nonlinear one. We assume 
that sign(O) = 0. There are three qualitatively different cases: A > 0, A < 0, 
and A = when we have the canonical model considered in [11,13]. 

2.2. Positive values of X 

Let us start from the case A > 0. It is convenient to set A = p 2 . The 
equation of motion for this case takes the form 

(^-^)0 + sign(0)-pV = O. (4) 

The potential has one local minimum at <f> = and two local maxima at 




Fig. 1. The potential V{4>) (A > 0) and its limit for A = 0. 

(f) = ±4j, see Fig.l. It is not differentiable at its minimum (right-hand side 
and left-hand side derivatives are not equal at this point). Equation (4) can 
be obtained from equation of motion which describes a small perturbation 
around the ground state in the model considered in [12]. Physical values of 
perturbation are given by \(p\. 

2.3. Negative values of X 

It turns out that the model for A < has also a physical meaning. Phys- 
ical system which is related to this model can be obtained from the system 
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considered in [11] by adding new springs that link every ball with the floor, 
see Fig. 2. In this system every ball can move only in vertical direction. There 




Fig. 2. The system of bouncing balls. 



is a rigid floor at <fr = and every ball bounces elastically from the floor. 
After taking few standard steps (continuous limit, folding transformation) 
we get the system which dynamics is described by the equation 

(df-^)</> + sign(0) + aV = O, (5) 

where — oo < <f> < oo whereas physical position of the balls are given by \4>\. 
As above, we have set A = —a 2 . The potential Fig. 3 has one minimum and 
no local maxima. 




Fig. 3. The potential V(<f>) (A < 0) and its limit for A = 0. 
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2.4- On regularized potentials 

In our model at 4> = first derivative of the potential does not exist at all. 
It is possible to replace the sharp potential by a regularized potential where 
instead of \d>\ we use e.g. \J § 2 + e 2 or e In cosh |, what gives well defined first 
derivative at <f> = 0. We are not interested in regularized potential in this 
paper because physical systems that we consider here give rather sharp than 
regularized potential. Another reason is that for the regularized potential 
only solutions such that \d>\ 3> e can survive the limit e — > 0. Finding any 
such solution of the equation (3) with the regularized potential seems to be 
a very difficult task. 



3. Symmetric shock waves inside the light cone 

3.1. The Ansatz 

Equation (3) includes term with derivatives in the form of d'Alembert 
operator so it can be reduced to an ordinary differential equation by assump- 
tion 4>{x,t) = W(z), where z = (x 2 — t 2 )/<i. 

There are two qualitatively different cases (z < and z > 0) for the 
canonical model (A = 0). When z < the pieces of solution can be combined 
together in one solution. In opposite case whole solution can be either posi- 
tive or negative and it is not limited respectively from above or from below. 
It is physically reasonable to get rid of unstable solutions. It can be simply 
done by modification our Ansatz to the following one 4>(x,t) = 0(— z)W(z), 
where @(—z) is well known Heaviside step function. This modification in- 
troduces discontinuities at x = ±t. 

In the modified model, especially for A > there are solutions for z > 
that are limited both from above and from below so this time we can not 
simply get rid of them. In further part of our work we analyze solutions for 
z < and z > separately because the solution inside the light cone is com- 
pletely independent of the solution outside it. It is executed by assumptions 
<f>(x, t) = @{-z)W(z) and <f>(x, t) = @{z)W(z). 

Apart from possibility of reduction equation (3) to an ordinary differ- 
ential equation, an important question is whether discontinuities in our 
model can move with velocity v ^ 1 or not. To answer it let us consider 
(f)(x,t) = Q(—z)W(z), where z = pq, p = (x — vt)/2, q = (x + vt)/2 and 
v 0. Equation (3) takes the form 



G(-z) 



A(p, q)W" -\{v 2 + l)(zW" + W) + sign(VF) - XW 



+ 



-25(z)A(p, q)W'(z) - 5'(z)A(p, q)W(z) = 0, 
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where 

A(p,q) = ^(v 2 -l)(p 2 + q 2 ), 

and ' = 4z. The terms proportional to @(—z), S(z) and S'(z) have to vanish 
independently. At z = and v ^ 1, W(0) = and VF'(O) = because 
A(p, q) 7^ 0. It means that (f>(x, t) can not be discontinuous function unless 
v = 1. We can see that velocity t> = 1 is distinguished by the model and it 
is the only admissible velocity with which discontinuities can move. 



3.2. Equations of motion 
Let us consider the Ansatz 

<f>(x,t) = Q(-z)W(z), where z = ^{x 2 -t 2 ). (6) 

Applying Ansatz (6) to equation (3) we get the following differential equation 
zW" + W - sign(VF) + XW = 0. (7) 
Let us introduce a new variable y which is related to z in the following way 




(8) 



Consequently, equation (7) acquires more familiar form 

G" + ^G'-p 2 G = -sign(G), A>0 (9) 

or 

F" + ^-F' + a 2 F = -sign(.F), A < 0. (10) 

Equations (9) and (10) are Bessel equations with the signum nonlinearity. 
We have denoted W(— \y 2 ) = G(y) for positive values of A and W(— \y 2 ) = 
F(y) for opposite case. 



3.3. Solutions for A > 

The term sign(G) is constant (equal ±1) on the intervals where sign of 
G(y) is constant. Equation (9) has the solutions: 

G + (y) = ^ - ah(py) - (3K (py) for G(y) > 0, 

G-(y) = ~ + alo(py) + (3K ( P y) for G(y) < 0, 
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where a, (3, a, (3 are arbitrary constants. We will find these coefficients 
matching pieces of solution so as to obtain a solution that is valid on full 
available range of y. Physically, only real solutions are interesting so only 
y > is considered. It is convenient to introduce the coefficients and (5k 
and combine all the solutions in one formula 



where k = 0, 1, 2, We assume that Go(y), that is G + type, starts from 

y = 0. It does not cause loosing of generality because physically relevant 
quantity is \<f>\, thus case where G-(y) plays role of solution that start from 
y = does not have to be extra discussed. For Go(y) the coefficient fio has 
to vanish to ensure regularity at y = 0. The coefficient «o can be expressed 
with help of Gq(0), and then Go(y) takes the form 



For fixed p there are three qualitatively different cases depending on Gq(0). 




(11) 




(12) 




4: 



X 



Fig. 4. The unstable symmetric shock wave at t=2. 




To 



Fig. 5. The unstable symmetric shock wave at t=5. 




Fig. 6. The unstable symmetric shock wave at t=12. 
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For Go(0) > \ the solution Go(y) — > oo for y — > oo (unstable solution). 
In this case solution Go(y) cover whole range y > 0. The shape of the shock 
wave for p = 0.5 and Gq(0) = 6.0 at different times is shown in Figs. 4-6. 

In the case Go(0) = \ we obtain the shock wave 

<f>(x,t) = -^9(t 2 -x 2 ). 

For this wave, values of the field behind the wave front are constant and 
equal l/p 2 . 

Last case Go(0) < jp is more complicated but it is much more interesting. 
The solution Go(y) holds only on the interval < y < 6o, where Go(y) > 0. 
For fixed Go(0), the first zero of G{y) i.e. bo is determined by solution of 
following equation 

h{pbo) 



i-Go(oy 



Unfortunately, it can be solved only numerically. It is clear that &o depends 
on p. For given bo first piece of solution G(y) takes the form 



Go(y) = j 



1 h(py) 



h{pbo) 



(13) 



We are interested in a solution for all nonnegative values of y. Having pieces 
of solution (11) and matching conditions 

G fc (6 fc _!) = 0, G" fe (& fe -i) = GLife-i), (14) 

which are implied by equation (9) we can calculate coefficients and 
First matching condition allows to eliminate coefficients Solution (11) 
takes the form 



G k (y) = (-1) 



i ( l K (py) 



P 2 \ Koipbk^ 

j(h \ ( J o(^) K o(py) \ 



(15) 



The zeros b^ for A; = 1, 2, . . . come from the equation Gk(bk) = 0. Like for 
k = we can get them numerically. The second condition in (14) gives 
coefficients a^: 

1 JC(b ,b ) - p 2 a h(pb ) 

(X\ = 

p 2 I (pb ) l(b , b ) + K.(b , b ) 
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where 



and k = 2, 3, . 



1 



Oik 



fc-i, 



A- p 2 I (pb k _ 2 )B 



A = 



B = 



/C(6 fe _i,6jk_ 2 ) +K,(b 



k-l,bk-l 



l(b k -i, b k -i) + /C(6 fc _i, & fc _i) 

^(6fc-l, &fc-2) + £(&fc-l, &fe-2) 



We have introduced the special notation 



h{pyY 



JC(x,y) = 



Kijpx) 




Fig. 7. The solutions G(y) for fixed Go(0) = 0.5 and different values of p. The curve 
^4 correspond to p = 0.1, B to p = 0.5 and curve C is the solution for p = 1.4. For 
fixed Go(0) the parameter p belongs to the interval < p < l/y / Go(0). 

The zeros of G(y) depend on p - they are larger for larger values of p, Fig. 7. 
For fixed Go(0) and p — > Go(0) -1 / 2 first zero bo — > oo and, of course, bk — > oo 
(bo < b\ < 62 < • • •)• All the zeros bk are larger than their counterparts in 
the canonical model. It means that zeros Xk = ± (t 2 — ^ 2 in the modified 
model run faster than x c k = ± (i 2 — a 2 ,) 1 '' 2 in the canonical model. A pair of 

zeros Xk appears at t = bk and moves with velocities v k = ± (1 — b\/t 2 ) 

In Fig. 8-10 we present three snapshots which show the symmetric shock 

wave for Gq(0) < ^j. 



5.^. Solutions for A < 

This section is devoted to presentation of solutions correspond to equa- 
tion (10). Many steps are the same so we sometimes skip the comments. 
Let us start from the solution of (10) given in the form 

F k (y) = (-if (pkJoW + VkY (ay) - j^j , (16) 
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X 



Fig. 8. 



-2-1 12 

The symmetric shock wave for p = 0.5, Go(0) = 0.5 at t=0.5. 





\ 0.4- 


l*l J 






\ 0.2: 




X 



Fig. 9. The symmetric shock wave for p = 0.5, Gq(0) = 0.5 at t=1.6. 





0.4: 


Itl j 






a 0.2- 




X 


-4 U 4 



Fig. 10. The symmetric shock wave for p = 0.5, Go(0) = 0.5 at t=4.7. 



where Jo and Yq are Bessel functions. They take real values for y > 0. F^(y) 
are positive for fe = 0, 2,4, . . . and negative for k = 1, 3, 5, . . .. We assume 
that Fo(y) starts from y = 0. This time, for given a there is no qualitatively 
change in behaviour of solution for different -Fb(O). Fo(y) can be expressed 
in the form 



F (y) = (f (0) + 4^ J (ay) - -1. 



(17) 



The first zero Co is calculated from the equation Fq(cq) = 0. We can rewrite 
F (y) using c 



Fo(y) = ± 



Jo(crco) 



1 



(18) 



In order to have solution for whole range y > we have to calculate /ik and 
ffc. We use the matching conditions 



^(Cfc-l) = 0, i£(c fc _i) = F^^Cfc-i). 



(19) 



For fe = 1,2, 



^(y) = (-i) 



fc+i 



^o(*l/) 



2 ^ YoK/i-l) 
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- fi k J (ac k -i) 



Yo(ay) 
Y (ac k -i, 



(20) 



The coefficients v k have been eliminated by using first condition in (19). 
The zeros c k fulfil the equation F k (c k ) = 0. Second condition in (19) gives 
coefficients \i k 



Mi = ~ 



1 y(cQ,C ) + H 2 wJl{pC 

a 2 J (ac ) J(cq, c ) - y{c , c ) 
1 



where k = 2, 3, . . . and 



a 2 Jo(ack-i) 



A + a 2 J (ac k ^ 2 )B 



A = 



B = 



yfa-i, c fc _ 2 ) + y(c fc -i, Cfe-i) 
j{c k -i-,c k -i) -y{c k -i,c k -i) 

J{c k -l,C k - 2 ) -y(Ck-l,C k -2) 



J{c k -1, Cfc_l) - y(Cfc_l, Cfc_l) ' 

By analogy, the functions and ^ are defined by formulas 



Ji (era;) 



y{x,y) = 



*b(<ri/)' 



Having /Xfc we can find consecutive zero as the solution of the equation 
F k (c k ) = 0. Fig. 11 presents solutions -F(y) for i*b(0) that is set and different 
values of a. The zeros c k are smaller than in the canonical model so zeros 



0.4' 



■ 










■ 


c\ b\a 

— i ■ \ i ■ — -^^W- — i — 








i ■ — 




/ 2 £ 







-0.2 

Fig. 11. The solutions F(y) for fixed Fq(0) = 0.5 and different values of a. The 
curve A correspond to a = 0.1, B to u = 0.5 and curve C to a = 2.0. 

Xfc = ±(1 — c\/t 2 ) 1 / 2 move slower than zeros x k . They are getting smaller for 
a being increased. The shape of \4>\ qualitatively resembles that presented 
in Fig.8-10. 
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3.5. Divergence of sequences bk and Ck 

An interesting problem is if the sequence of bk (or c&) is divergent or not. 
Unfortunately, we can not show a proof of divergence these sequences (only 
a numerical evidence) because explicit expressions for bk (ck) are not known 
- we have got them only as a result of numerical computations. To see it, 
we introduce 



b n = b xix 2 



where 



Xk 



bk 
bk-i 



Let us say that for n — > oo the sequence of b n is divergent b n — > oo. It also 
means that the sequence of In b n is divergent and sum of these logarithms 



i=i 



as well. We apply Kummer criterion so as to check if the series is divergent 
[16]. In our computation the comparative sequence q n = nlnn has been 
used. It has been checked up to n = 700 that 



K-n — Qr 



b n 
b n +l 



Qn+1 



are negative and monotonically decrease, see Fig. 12. Kummer criterion says 




] 100 200 300' ' ' 400' ' ' 500' ' ' 600' 

Fig. 12. Behaviour of tC n for sequence b n . 



that series is divergent if for all n > N, values of JC n < 0, where N is a fixed 
number. It suggests that solution G(y) covers whole y > 0. It has been 
checked that also the sequence of c n is probably divergent. 



3.6. The correspondence between the modified and the canonical model 

The potential V(<f>) = \<f>\ is the limit case of the modified one for A — > 0. 
We are interested in limit A — > for the solutions in the modified model. It 
is not a priori clear that these limit solutions have to be solutions that are 
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known from the canonical model. We are able to do an analytical computa- 
tion for first two pieces of solution i.e. Go and G\ (or Fq and Fi). Let us 
remind that two first pieces of solution W c (z) for the canonical model have 
the form 

Wl 1 (z) = z + a , W$(z) = - ( z + a + 2a In ^ V (21) 

V «o/ 

where — ao is the first zero (for variable z). For more detail see [6]. Let us 
denote WZ^z) = G (y(z)) and W + (z) = Gi(y(z)) (by analogy WZ^z) = 
Fo(y(z)) and W ~(z) = Fi(y(z))). We also rename p + = p and p_ = a. We 
are interested in comparison of the canonical and the modified solution for 
p± — > 0. In this limit the condition Wl^O) = WZ~i{0) can be replaced by 
ao = clq , where = bo/4 and Oq = cjj/4. The solution and can 
be expanded in the Taylor series what gives 

W^(z) =z + a TCp 2 ± + O(pi), (22) 



where 



W^(z) = - (z + a + 2a ln^) T Dp 2 ± + O(pi), (23) 



z 2 



^ z ° 2 

C = — + a 2;+ -a , 



D = (ag - 2a ^) In + 3a .z + -rOo- 

ao 4 4 

The leading terms in (22) and (23) do not depend on p± so they are limits 
of these solutions for p± — > 0. The most important thing is that these limits 
are exactly equal to the first two solutions WH 1 and W§ in the canonical 
model. The correspondence between the other solutions can be checked 
numerically. We can see that the smaller values of k we take, the better 
correspondence is. 

4. Symmetric shock waves outside the light cone 

4-1. Equations of motion 

The solutions outside the light cone can be obtained with help of the 
Ansatz 

(j)(x,t) = 0(z)W(z), where z = ^(x 2 -t 2 ). (24) 
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If we now introduce z = \y 2 we get 



g" + -g' + p 2 g = sign (5), 



(25) 



r + -/'-^ 2 / = sign(/), 

y 

where g{y) = W(z(y)) for A > and f(y) = W(z{y)) for A < 0. 



(26) 



4.2. Case A > 

It is convenient, for our further analysis of solutions, to associate the 
potential U(g) with equation (25). We can easily see that U(g) = \p 2 g 2 — \g\, 
see Fig. 13. 





U(g) 
















Fig. 13. The potential U(g). 
The solution of (25) takes the form 

9k(y) = (-!) fe (^2 - PkMpy) - VkY {py)j , 



(27) 



where g^ > for k = 0, 2, 4, . . . and gt < otherwise. As above, we have to 
set uq = so as to have go(y) regular at y = 0. If go(0) is given then 



go(y) 



l 



50 (0) J (py). 



(28) 



In Fig. 14 we present a few curves (28) for different values of <?o(0). There 
are several qualitatively different cases. 

Case <7o(0) = corresponds to solution that starts from the point where 
U(g) has its local maximum. We have chosen positive solution but for 
<7o(0) = 0, negative solution is possible as well. We can get these solutions 
by replacing (— l) k in (27) with (— l) h+1 and assuming that go(y) starts from 
y = 0. We will not discuss this situation separately because \<j>\ is physical 
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te- 

14- 

12: 

8^ 




25 



Fig. 14. The curves go(y) for different go(0) and p — 0.5. They are solution for full 
range y > except those which start from go(0) = 14 and go(0) — 16. 

quantity. It is worth mentioning that this solution is not exactly shock type 
because there is not sharp front wave. 



the solutions go (y) are valid for all y > 0. j\ is the first zero of J±, (j® = 0). 
Approximately, p 2 g$ r ' it = 3.482872. This case contains the constant solution 
9o(y) = ji- 

For critical value of <?o(0) = g^ %t the first zero of g(y) appears. The 
solution can take either of forms: the solution 



for all y > or g (y) given by (29) for < y < j\/p and gi(y) = -g (y) for 

y > Ji /p- 

If go(y) is a little bit bigger than g^" 11 the solution is made up of go(y) 
and gi(y) but this time j\ in (29) is replaced with pco and g±(y) contains also 
function Yo(py), see Fig. 15. This picture is valid until go(0) reach consecu- 
tive critical value fffj"* 2 (unfortunately, we are not able to give appropriate 
analytical formula - only numerical value of g crit2 is available) . 

Therefore, for go (0) = <?o™* 2 situation resembles that for <?o(0) = go %t 
but this time two zeros Co and c\ already exist. 

For values of go (0) a little bit more than <?o™* 2 solution looks like it was 
shown in Fig. 16. For greater values of <?o(0) more zeros appear, see Fig. 17. 
An important issue is that for finite <?o(0), maximal value of k is always a 
finite number. For these solutions, asymptotic values of g(y) for y — > oo 



For < g (0) < g§ 



, where 





(29) 
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Fig. 15. The solution g(y) that has only one zero. 




Fig. 16. The solution g{y) with two zeros. 

















y 




j ' 20V / 4 







Fig. 17. The solution g(y) with six zeros. 



i.e. ±l//0 2 correspond to the minima of the potential U(g). We have not 
discussed yet the formulas for coefficients pk and Vk in (27). Fortunately, 
they have the same form as analogical coefficients in equation (16) - only a 
has to be replaced with p. 



4.3. Case A < 

The potential U(f) for equation (26) takes the form U(f) = ~^cr 2 f 2 -\f\ 
what suggests unstable behaviour of f(y). The solutions of (26) can be 
written down in the form 

f+(y) = al (py) + 0K o (py) - for f{y) > 0, 

f-(y) = -al (py) - (3K {py) + -L for f(y) < 0. 
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Let us consider f(y) > 0. By analogy to our previous analysis we will denote 
it as fo(y). For given /o(0) it takes the form 

fo(y) = (^ + /o(0))/ (<ry)- jjj, (30) 

where (3 = 0. There are five such solutions in Fig. 18. The solution fo(y) 
covers whole range y > 0. We have an analogical situation for a negative 
solution. The border solution (/o(0) = 0) is, of course, a non-shock type. 




3 1 2 3 4 5 6 7 



Fig. 18. The solutions fo(y) for different (positive) values of /o(0). 



5. Summary 

We have presented the shock wave solutions in the model with the po- 
tential V(4>) = \4>\ — \X4> 2 , where A is a real constant. The square term plays 
the role of perturbation of the potential V{4>) = \<j>\. The potential \<f>\ plays 
a special role because the corresponding equation of motion has the scaling 
symmetry. The square term is the simplest one that breaks this symmetry. 
It has been shown that both cases with nonzero values of A (A > and 
A < 0) have physical applications - appropriate potentials appear in the 
system of coupled bouncing pendulums or bouncing balls. 

Two kinds of waves have been found. The first one is exactly zero outside 
the light cone and has two wavefronts (the field is non-continuous at them) 
exactly at the surface of the light cone. Depending on Go(0) the solution 
inside the light cone is unstable, constant, or has isolated zeros. We have 
found that these zeros run faster (A > 0) or slower (A < 0) than their 
counterpart in the canonical model (A = 0). It has been also argued (by 
showing the numerical evidence) that zeros of solution, that depends on 
variable y, form probably divergent sequence. Moreover, we have shown 
(analytically in the case of first two pieces of solution) that shock waves 
inside the light cone for the modified model in the limit A — > reduce to 
solutions known from the canonical model. 

The second type of solution has also wavefronts at x = ±t. This solution 
takes zero values inside the light cone. There is family of solutions that 
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asymptotically (y — > oo) reach values ±l/p 2 . At the points 5 = ±l/p 2 the 
potential t/(g) (Fig. 13) has the local minima (or V{4>) has its local maxima, 
Fig.l). Among solutions that belong to this family, there are solutions that 
have no zeros, have only one zero, exactly two zeros etc. In contrast to 
the shock waves inside the light cone that have infinite number of zeros, 
solutions considered here have always finite number of zeros. There is also 
non-shock type solution - it has only one zero at the surface of light cone. 
It can be regarded as the border case of solutions belonging to this family. 
We have also found another family of solutions that have no zeros and grow 
to infinity for y — > 0. 
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